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Abstract — We consider a multiple-input multiple-output 
(MIMO) multiple access channel (MAC), where the channel 
between each transmitter and the receiver is modeled by the 
doubly-scattering channel model. Based on novel techniques from 
random matrix theory, we derive deterministic approximations 
of the mutual information, the signal-to-noise-plus-interference- 
ratio (SINR) at the output of the minimum-mean-square-error 
(MMSE) detector and the sum-rate with MMSE detection, which 
are almost surely tight in the large system limit. Moreover, we 
derive the asymptotically optimal transmit covariance matrices. 
Our simulation results show that the asymptotic analysis provides 
very close approximations for realistic system dimensions. 

I. Introduction 

Most works on wireless multiple-input multiple-output 
(MIMO) systems share the underlying assumption of a rich 
scattering environment and, thus, Rayleigh or Rician fading 
channel matrices with full rank. However, several measure- 
ments of outdoor MIMO channels have shown that this as- 
sumption fails to hold in certain scenarios, where low-rank 
channels are observed despite low antenna correlation at the 
transmitter and receiver (see e.g. JT|, 0). Motivated by these 
observations, a generalized fading MIMO channel model, the 
so-called "doubly-scattering model" (3), was proposed and 
has since then attracted significant research interest. A special 
case of the doubly-scattering model is the keyhole channel 
PI . J5] which exhibits null correlation between the entries 
of the channel matrix but only a single degree of freedom. 
The existence of such channels in reality was confirmed by 
measurements in |5|. 

Several theoretical works have studied the doubly-scattering 
model so far. The authors of [6| derive capacity upper-bounds 
for the general model and a closed-form expression for the 
keyhole channel. An asymptotic study of the outage capacity 
of the multi-keyhole channel was presented in Q. The diver- 
sity order of the doubly-scattering model was considered in 
|H) and it was shown that a MIMO system with t transmit 
antennas, r receive antennas and s scatterers achieves the 
diversity of order trsj max(£, r,s). A closed-from expression 
of the diversity-multiplexing trade-off (DMT) was derived 
in (9|. Beamforming along the strongest eigenmode over 
Rayleigh product MIMO channels, i.e., the doubly-scattering 
model without any form of correlation, was considered in 1 10 1. 
Here, the authors derive exact expressions of the cumulative 
distribution function (cdf) and the probability density function 
(pdf) of the largest eigenvalue of the Gramian of the channel 
matrix and compute closed-form results for the ergodic capac- 
ity, outage probability and signal-to-noise-plus-interference- 



ratio (SINR) distribution. In a later paper |QT|, the MIMO 
multiple access channel (MAC) with doubly-scattering fading 
is analyzed. The authors obtain closed-form upper-bounds on 
the sum-capacity and prove that the transmitters should send 
their signals along the eigenvectors of the transmit correlation 
matrices in order to maximize capacity. 

Despite the significant interest in the doubly-scattering 
channel model, little work has been done to study its asymp- 
totic performance when the channel dimensions grow large. 
We are only aware of 0, in which a model without transmit 
and receive correlation is studied relying on tools from free 
probability theory. Implicit expressions of the asymptotic mu- 
tual information and the SINR of the minimum-mean-square- 
error (MMSE) detector are found. 

In this paper, we consider a MIMO MAC with double- 
scattering fading in its most general form and derive deter- 
ministic approximations of the (ergodic) mutual information, 
the (ergodic) sum-rate with MMSE detection and the SINR 
at the output of the MMSE detector. The approximations 
become almost surely exact as the dimensions of all channel 
matrices grow large and can be easily numerically computed 
with negligible computing complexity. In addition, we provide 
the asymptotically capacity maximizing transmit covariance 
matrices and present an iterative water-filling algorithm for 
their computation. Our numerical results suggest that the 
asymptotic approximations are already very tight for channel 
dimensions with as little as four transmit and receive antennas 
and are therefore of clear practical value. 

The key idea behind the proofs in this paper is that 
the doubly-scattering channel model can be interpreted as a 
Kronecker channel lfl2ll with a random receive correlation 
matrix, which itself is modeled by the Kronecker model. This 
observation allows us to build upon [12] which provides an 
asymptotic analysis of the Kronecker channel model with 
deterministic correlation matrices. We then extend this work 
by allowing the correlation matrices to be random. The results 
in this paper are obtained through advanced tools from random 
matrix theory (inspired by fl3l . fBl . see also the textbook 
lTT3Tl for a comprehensive introduction and a contemporary 
overview of recent research results) and are hence not only a 
novel contribution to the field of wireless communications but 
also to the field of large random matrix theory. We also believe 
that the developed techniques can be successfully applied 
to the study of even more involved channel models, such 
as channels with line-of-sight (LOS) components or MIMO 
product channels with an arbitrary number of matrices. 



II. System model 

Consider a discrete-time MIMO channel from K transmit- 
ters, equipped with n k (k = 1, . . . , K) antennas, respectively, 
to a receiver with N antennas. The channel output vector 

y e 



at a given time reads 

K 

y = H fc x fc 



n 



(1) 



fc=l 



where U k e C Nxn " and x fc = [x kA , . . . , x k ^ nk } T ~ 
CA/~(0, Qfe), Qfc € C nfcX " fe , are the channel matrix and 
the transmit vector associated with the kth transmitter, and 
n ~ CA/"(0, pljv) is a noise vector. The channel matrices Hfc 
are modeled by the double-scattering model Q 

Hfc = i H.,-W;. A S; W 2 .fIV ( 2 ) 



^A^ 

iATxJV 



where R fc e C AJV , S fe e C***** and T fc e 
are deterministic correlation matrices, while Wi fe £ 
and W2,t € ^Nfcxrifc ^ ^dependent standard complex 
Gaussian matrices. Since the distributions of Wi^ and Wa,*, 
are unitarily invariant we can assume without loss of generality 
Sfc = diag(sfc j,..., s/j TVfc) to be diagonal matrices. Denote 
In(p) the instantaneous normalized mutual information of the 
channel ([T| in nats/s/Hz, defined as 



life xn fc 
,NxN k 



In(p) = ^logdet 



III. Main results 



(3) 



The notation N 
N k , n k grow infinitely large, satisfying 



00 denotes in the sequel that N and all 

lim inf ^ < 

lim sup ^ < 00, < lim inf ^ < lim sup t& < 00. 
These conditions ensure that all matrix dimensions grow at 
a similar speed. Additionally we need the following technical 
assumptions: 

A 1: For all k, limsup||Rfc|| < 00, limsup||Sfc|| < 00 
and limsup||TfcQfc|| < 00, where ||-|| is the spectral norm. 

Our first theorem introduces a set of 3K implicit equations 
which uniquely determines the quantites (g k ,9k, <5fc) (1 < k < 
K). These quantities will be used in the sequel to provide 
deterministic approximations of In(p) which become almost 
surely arbitrarily tight as N — > 00. 

Theorem 1 (Fundamental equations): The following sys- 
tem of 3K implicit equations in g k , g k and 5 k (1 < k < K): 



9k 



1 

nk 
1 

n k 



-(rT ; -Q,T| (g k T~;Q,Tl 



Sk,jSk 



S k = 



N, 



trR t 



I- 9kSk,jS k 

' K 

n k 9k9k 



(4) 



Rfc + pl 



N 



has a unique solution satisfying g k ,g k ,S k > for all k and 
p> 0. 



Remark 3.1: One can also prove that g k , g k and 5 k can 
be computed by a classical fixed-point algorithm which itera- 
tively computes Q, starting from some arbitrary initialization 
§k,gk,Sk > 0. This algorithm generally converges in a few 
iterations (depending on the system size) and does not pose 
any computational challenge. 

The next theorem provides a deterministic, asymptotically 
tight approximation of the (ergodic) mutual information based 
on the quantites (g k ,g k ,S k ) as provided by Theorem[T] 



Theorem 2 (Mutual information): 
(i) I N (p)-T N (p) 







(it) EI N (p)-I N (p) 



-> 



where 

In(p) 



N 



logdet I 



JV 



K 

E 



nk gkgk 



P-N k S k 



Ra 



1 K 

— ^2 [logdet (I Nk + g k S k S k ) 



fc=i 



+ logdet (I nk + g k T* Q fc T fc 2 j - 2n k g k g k 
and g k ,g k ,S k are the unique positive solutions to Q. 



The following result allows us to compute the asymptoti- 
cally optimal precoding matrices Qfc which maximize In(p) 
under individual transmit power constraints. 

Theorem 3 (Optimal power allocation): The solution to 
the following optimization problem: 



(Qi, ,Q 



K 



arg max I N (p) 
Qii— iQk 

s.t. — trQ fc < P k Vfc 

n k 

is given as Q* = U fe P£U£, where U fc <= C" feXnfc 
is defined by the spectral decomposition of Tfc = 
U fe diag(i M ,...,ifc,„JU£ and P* = diag(p fe l , „ J is 

given by the water-filling solution: 



Pk, 



Pk 



9pK 



where p k is chosen to satisfy — trP^ 
by Theorem [l] for Qfc = Q*. 



(5) 

Pk and g* k is given 



Remark 3.2: The optimal power allocation matrices P* k can 
be calculated by the iterative water-filling Algorithm 1 (see 
lfl2l Remark 2] and lfl3l Remark 5] for a discussion of the 
convergence of this algorithm). 

The last two results of this correspondence provide de- 
terministic approximations of the SINR at the output of the 
MMSE detector and the sum-rate with MMSE detection. 



Algorithm 1 Iterative water-filling algorithm 



4: 



Let e > 0, n = and p* k °j = P k for all k, j. 
repeat 

For all k, compute g* k n according to Theorem [T] with 
matrices Q k = U fc diag (AjJ U£. 

For all k,j, calculate p* k '™ +1 = (fJ> k 

n k such that ^E"ii^;" +1 = p ^ 



with 



5: n = n + 1 
6: until maxtj 



-*,n— 1 



< e 



Theorem 4 (SINR of the MMSE detector): Assume Q k = 
I nk and T k = diag(7 M , . . . ,t k<nk ) for all k and let 7 fej be 
the SINR at the output of the MMSE detector related to the 
transmit symbol x k j, given by 



K 



7*^ = h L E 11 "' - h *Xi + ?M h 



v»=l 



Then 



7fe,j - 7fc„ 



N-yoo 







where 7^ = t k jg k and ^ is by given by Theorem [T] 

Remark 3.3: The theorem is also valid under the more 
general assumptions T fe = U fc diag(f M , . . . , t k , n k )V k and 
Q fc = U fe diag(p fcj i, . . . ,Pk,n k )^k- We can then sim p!y de- 
fine the matrices T' k = di&g(t kj ip k ,i, . . . ,tk,n k Pk,n h ) and 
Q' fc = I]y k for which the theorem holds. 

Corollary 1 (Sum-rate with MMSE decoding): Under the 
same assumptions as in Theorem |4] let 



1 _ K . J lk . 
fe=li=l 



' 7fc, 



Then, 



K n h 



-> 



fc=lj=l 

(it) Ei?(p)-l^£log(l + t fcjm ,)-^0 

iv * — ' ^ — ' AT— foo 

k=U=l 

where g k is given by Theorem [T] 

A. 77ze Rayleigh product channel 

A special case of the double-scattering channel is the 
Rayleigh product MIMO channel [10] which does not exhibit 
any form of correlation between the transmit and receive 
antennas or the scatterers. For this model, the Theorems [T] 
[2] and |4] can be given in closed-from as shown in the next 
corollary. 



Corollary 2 (Rayleigh product channel): For all k, let 
N k — S, n k — N and assume = Ijv, S k = Is and 
Rfe = Iff. Then, 

7jv(p) = log 1 + --^-5 15+^-1 



iV 



and 



log 1H — fa - 1 

N 6 V 5 

X log (3) -2Jf(l-g) 
1-5 



7fc,j = 



9 



where g is the unique solution to 



9(1 



s 


1 


N 


~ K 


S 


1 


N ' 


~ K 



s 



(1 + p) 



5 



(6) 



such that 5 = (1 - g)/(g(g + S/N - 1)) > and = (1 - 

Note that similar expressions for the asymptotic mutual 
information and MMSE-SINR have been obtained in Q by 
means of free probability theory. However, their results require 
the numerical solution of a third order differential equation. 

IV. Numerical examples 
As a first example, we consider the "multi-keyhole channel", 



i.e., N k = 1, S k 



1, 



I nk , for 



N = n k = 4. The signal-to-noise-ratio (SNR) is denoted by 
SNR = 1/p. Fig. [T] depicts the normalized ergodic mutual 
information Eijv(p) and its asymptotic approximation In(p) 
versus SNR for a single (K = 1) and multiple (K = 3) 
"keyholes". Surprisingly, the match between both results is 
almost perfect although the channel dimensions are very small. 

As a second example, we consider a multiple access channel 
from K — 3 transmitters, assuming the double-scattering 
model in |3| . Under this model, the correlation matrices are 
given as R fe = G(<f> rik ,d r>k ,N k ), S k = G((/) Stk ,d Stk ,N k ) and 
T k = G(<f>t,k,dt t k,Nk), where G(<£,d,n) is defined as 



[G(<f> ) d ) n)] k>l = - E exp fi27rd(fc-Z)siiW - 



The values 4> t , k and </> rj fc determine the angular spread of the 
radiated and received signals, d tyk and d rk are the antenna 
spacings at the fcth transmitter and receiver in multiples of 
the signal wavelength, and N k can be seen as the number and 
d s>k as the spacing of the scatterers. For simplicity, we assume 
N = 4, P k = l/n k , N k = 11, n k = 3, d t , k = d r , k = 0.25 
and d s>k — 50 for all k. We further assume <f> r>k — (j> t ,k for 
all k, with (f> r ^k € {tt/4, 7r/2, tt} and Si fe = 7r/8. Fig. [2] 
shows E/jv(p) and In(p) with uniform and optimal power 
allocation versus SNR. Again, our asymptotic results yield 
very tight approximations for even small system dimensions. 
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Fig. 1. Ergodic mutual information EJjv(p) of the multi-keyhole channel 
and its asymptotic approximation In(p) versus SNR. 
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Fig. 2. Ergodic mutual information Eijv(p) an d sum-rate Eijjv(p) of 
the multiple access channel and their asymptotic approximations /jv(p) and 
Rn(p) versus SNR. 



optimal transmit covariance matrices. Our numerical results 
show that the asymptotic analysis provides very close approx- 
imations for very small system dimensions with as little as four 
transmit and received antennas. We believe that the techniques 
used in this paper could be succesfully applied to the study of 
even more involved channel models. 

VI. Appendix 

Sketch of proof of Theorem [7} The key idea is that the 
doubly scattering model can be considered as the Kronecker 
channel model [12| with random correlation matrices. For the 
Kronecker model, the matrices are given as 

1 



H 



k — 



:Z fc Wo fe T 



(7) 



where Z k e C Nx k is a deterministic matrix and \V2.fc 
and Tfe are as defined in The fundamental equations (cf. 
Theorem [T| for this model are given by lfl2l Corollary 1] 

1 1 1/1 1 \ -1 



ejt = — trZ fe Z£ ( Ve.Z.Z," + pl N ) 

Assume now Zj. to be random and modeled as 

1 



Z fc = 



(8) 



(9) 



Notice first that the expressions of the quantities e k are un- 
affected by this assumption. Second, e k have become random 
quantities and it is our goal to find deterministic approxima- 
tions g k of e k , such that e k — g k as N — > 00. Following 
similar steps as in the proof of IPT31 Theorem 4], one can now 
show that 



max |efc 

k 



max \ ek 

k 



9k\ 



9k\ 



N- 



-> 



-> 



(10) 



where g k and g k satisfy Q. The proof of uniqueness of such 
solutions relies on arguments of so called standard functions 
lfT6ll and follows similar steps as in |[T3l Proof of Theorem 3] 
or (Hi Proof of Theorem 1], ■ 



For comparison, we also provide the sum-rate with MMSE 
detection EiJjy(p) and its deterministic approximation R^. 
We observe a good fit between both results at low SNR 
values, but a slight mismatch for higher values. This is due 
to a slower convergence of the SINR j k j to its deterministic 
approximation 7^ . 

V. Conclusion 

We have studied a MIMO MAC with doubly-scattering 
fading channels. Under the assumption that the dimensions 
of all channel matrices grow infinitely large, we have derived 
almost surely tight deterministic approximations of the mutual 
information, the SINR of the MMSE detector and the sum- 
rate with MMSE detection. In addition, we have provided an 
iterative water-filling algorithm to compute the asymptotically 



Sketch of proof of Theorem |2j We rely again on the 
observation that the doubly scattering model can be considered 
as the Kronecker channel model [12] with random correlation 
matrices Z k (cf. (|7]i and Q). A deterministic equivalent 
lff on (p) of the mutual information for the Kronecker model 
(|7| was provided in [12, Theorem 2]. Here, J^ on (p) is defined 
as a function of the quantities e k and e k (given as the unique 
solutions to ([8]l) which reads: 



fKron 

N 



(P) 



N 



log det 




K 1 

Y- 

k=l 



log det (l nk +e fe T*Q fc T* 



1 K 

-Y 

k=l 



nk^ke-k- 



(11) 



k 



Due to the almost sure convergence of e k — g k and 
efc — 9k — > established in Theorem [T] we can simply 
replace e k and e k by <?fc and g^, respectively. It remains 
now to find a deterministic equivalent of the first term 
~ logdet + ~ Sil-Li CfcZfcZ^V which is random since 
the matrices are random. However, this term is nothing 
but the mutual informatipn of another Kronecker channel with 
matrices H 



-R^Wi^Sj . Hence we can apply again 
[12 Theorem 2] to find its deterministic equivalent. Combin- 
ing both results yields In(p) and terminates the proof of (i). 
Denote fl the probability space engendering the sequences 
{Wi t i, . . . , Wi jfjWj i, . . . , W2,jf}. Then, on a sub-space 
of £1 of measure 1, we have by (i): In(p) — In{p) — > as 
N — > oo. Integrating this expression over £1 implies by the 
dominated convergence theorem ifTTl part (ii). ■ 

Sketch of proof of Theorem [if Similar to the proof of 
lfl2l Proposition 3], one can show that the covariance matrices 
Q* k should align to the eigenvectors of the transmit correlation 
matrices to maximize In(p), i-e-, Q* k = U^P^U^. Note 
that it was also proved in IfTTl Theorem 1] that these signaling 
directions are optimal to maximize EJjv(p). One can then 
further show that 



d Pl, 



9*k t k,j 



1 + 9Pk,jPt 



Wk,j 



and 



d 2 I N (p) 



< 0. Since In{p) is hence strictly concave in p* k ■ 

it follows from the KKT conditions lfl8l that P? is given by 
the water-filling solution |5) with power constraint ^-trQ£ = 

&n = Pk- ■ 

Sketch of proof of Theorem^ Assume that H^, are given 
by the Kronecker model (|7]i with deterministic matrices Z k . It 
follows from standard lemmas of random matrix theory (see 
e.g. Ifl3l Appendix C], that the following limit holds: 



t fc)J -^-trZ fe Z£ ^H,H,H+pI„^ 



0. 



Applying [12 Theorem 1] to the second term yields 



fcj — trZ fc Z£ (f>Z,Zr+pI„ fc ] -5^0 



where are given as the solutions to (|8). Notice from ([8]) that 
the second term is equal to t k je k . Assume now the matrices 
Z k random and given by Q. Thus, we have by ( fT0| ) 



Jk,j — tk,j9k 



N- 



-> 0. 



Sketch of proof of Corollary^ Part (z) follows directly 
from Theorem [4] and the continuous mapping theorem [19 
Theorem 2.3]. Part (ii) follows from the same arguments as 
in the proof of Theorem [3] (ii). ■ 



Sketch of roof of Corollary^ One can show by straight- 
forward but tedious calculations that the fundamental equa- 
tions Q can be reduced to a single implicit equation |6]) if 
T k = Ijv, S fe = I s and K k = I N for all k. Replacing g k 
and 5 k by g and 5 in the expressions of In(p) and jk,j, 
respectively, leads to the desired result. ■ 
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